s-wave superconductivity probed by measuring magnetic penetration depth and lower 
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The magnetic penetration depth A has been measured in MgCNi3 single crystals using both a high 
precision Tunnel Diode Oscillator technique (TDO) and Hall probe magnetization (HPM). In striking 
contrast to previous measurements in powders, <5A(T) deduced from TDO measurements increases 
exponentially at low temperature, clearly showing that the superconducting gap is fully open over the 
whole Fermi surface. An absolute value at zero temperature A(0) =230 nm is found from the lower 
critical field measured by HPM. We also discuss the observed difference of the superfiuid density 
deduced from both techniques. A possible explanation could be due to a systematic decrease of the 
critical temperature at the sample surface. 

PACS numbers: 74.25.Nf, 74.25.Op, 74.70.Dd 



The interplay between magnetism and superconduct- 
ivity is currently a subject of great interest. In the UGe2 
and URhGe uranium compounds, for instance, a long 
range ferromagnetic ordered phase coexists with the su- 
perconducting phase and a mechanism of spin fluctu- 
ations (SF) could be at the origin of the Cooper pair 
formation [l|, 0|. The recent discovery of high temper- 
ature superconductivity in oxypnictides also rapidly be- 
came the topic of a tremendous number of both experi- 
mental and theoretical works. The parent undoped LnO- 
FeAs (where Ln=La,Sm,...) compound is here close to it- 
inerant magnetism due to the presence of a high density 
of Fe d states at the Fermi level [U, leading to competing 
ferromagnetic and antiferromagnetic fluctuations. Sim- 
ilarly in the cubic (anti-)perovskite MgCNi3 compound 
[3], the presence of a strong Van Hove singularity in the 
density of Ni states slightly below the Fermi level also 
leads to strong ferromagnetic fluctuations These 
two systems have also a Fermi surface composed of both 
electron and hole pockets (3D sheets in MgCNi3 as com- 
pared to quasi-cylindrical sheets in oxypnictides). 

Despite these striking similarities in their electronic 
and magnetic properties, spin fluctuations lead to very 
different effects in those systems. On the one hand, 
ab-initio calculations rapidly showed that the electron- 
phonon coupling constant (X e - P h ~ 0.2) is far too low 
to account for the high critical temperatures observed 
in oxypnictides (up to ~ 55 K) and an unconventional 
mechanism mediated by the SF associated with a sign 
reversal of the (s-wave) order parameter between elec- 
trons and holes sheets of the Fermi surface has been pro- 
posed [101 . On the other hand, it has been suggested that 
the narrow van Hove singularity could be responsible for 
a nearly unstable phonon mode in MgCNi3 inducing a 



large, although reduced by SF, \ e ~ p h [ill . fl2| in agree- 
ment with experiments which yield an average electron- 
phonon coupling constant in the order of 1.7 jlll [13, fl3| . 
The interplay between electron-phonon coupling and SF 
is further emphasized in this system by the existence of 
a large isotopic effect [l5[ which has been suggested to 
be enhanced by the strong SF (l6j . 

In this context, the nature of the superconducting 
order parameter rapidly became a crucial issue. In 
MgCNi3, the experimental results still remain contro- 
versial : on the one hand, penetration depth measure- 
ments (in polycrystalline samples) showed a quadratic, 
i.e. non s-wave, temperature dependence suggesting a 
nodal order parameter [17| . whereas specific heat meas- 
urements clearly indicate that the superconducting gap 
(A) is full y op en, with a A/fcsT c ratio ranging from 1.9 
to 2.1 (l3l.ll8l. [l9j i.e. well above the BCS weak coupling 
1.76 value. 

In this paper, we present high precision magnetic pen- 
etration depth and lower critical field measurements per- 
formed in the same MgCNi3 single crystals. We show 
that A(T) clearly follows an exponential temperature de- 
pendence for T < T c /3 showing that the gap is fully open 
on the whole Fermi surface. A zero temperature Ao value 
of 230 nm, i.e. well above the London clean limit BCS 
value (~ 60 nm) has been deduced from first penetra- 
tion field measurements, clearly suggesting the presence 
of strong mass renormalization and/or impurity scatter- 
ing effects. Introducing this value into the TDO data 
however leads to a temperature dependence of the nor- 
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is different from the one directly deduced from the lower 
critical field (H c i oc Iti(k)/X 2 where k = A/£). Possible 
reasons for this discrepancy are discussed. 
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Single crystals were grown in a high pressure furnace 
as described elsewhere j20|. AC specific heat have been 
performed on several samples of the same batch plj . 
All the measured crystals show sharp superconducting 
transitions (AT C ~ 0.2 K) emphasizing the excellent bulk 
homogeneity of each crystal. We observed however a 
large dispersion of critical temperatures from sample to 
sample, between approximately 5.9 to 7.6 K, probably 
due to a slight Ni deficiency in the MgCNi3 structure 
(20j . Three single crystals with a thickness of 0.1mm 
but different shapes and critical temperatures have been 
selected. Sample ^A can be approximated by a disk 
with a diameter of 0.3mm. Samples #B and #C have a 
rectangular shape of 0.21x0.15 mm 2 and 0.24x0.36 mm 2 
respectively. Samples #A and #B both present a bulk 
T c of 6.9 K and exhibit exactly the same behavior by 
TDO and HPM, whereas sample #C has the highest T c 
at 7.6 K. 

The magnetic penetration depth has been measured 
with a high stability LC oscillator operating at 14 MHz, 
driven by a Tunnel Diode [22I |23| . The AC excitation 
field is below 1 pT and the DC earth magnetic field is 
screened by a demagnetized weak ferromagnet amorph- 
ous ribbon, ensured to work well below H c \. The sample 
stage, placed at the bottom of a home-made He 3 refri- 
gerator, is regulated between 0.5 K and 10 K, whereas the 
LC oscillator remains at fixed temperature. The super- 
conducting sample is glued with vacuum grease at the 
bottom of a sapphire cold finger, which can be extrac- 
ted in-situ (24|. The small filling factor of the excitation 
coil by the superconducting sample (~ 0.01%) ensures a 
small perturbation of the circuit and the frequency shift 
Sf divided by the one induced by the extraction of the 
superconducting sample, A/q, is then proportional to the 
imaginary part of the surface impedance and hence to the 
magnetic penetration depth |25|. As shown in the inset 
of Fig. 1, all samples present a sharp superconducting 
transition at a critical temperature T[ (defined by the 
onset of the frequency shift change) equal to 6.9 K (resp. 
7.6 K) for sample #A (resp. sample #C) in good agree- 
ment with the C p measurements. 

Fig.l displays the temperature dependence of the fre- 
quency shift, proportional to (5A(T), compared to the res- 
ults previously reported in powders (l7| . The amplitude 
of the shift is 10 times larger in the case of the powder 
reflecting the fact that the surface on which the super- 
currents are flowing is much larger in powders than in 
single crystals (for the same sample volume). It is im- 
portant to note that the temperature dependence of A 
is strikingly different in single crystals than in powder 
for which a T 2 power law has been reported below 1.8 K. 
Such a dependence has been interpreted as an evidence 
for unconventional superconductivity fl7| but our meas- 
urements do not support this scenario as a T 2 power law 
only very poorly describes the experimental data (see 
dashed line in Fig.l). 

A very good fit to our data is actually obtained assum- 
ing the low temperature approximation for clean type II 
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Figure 1: Low temperature TDO frequency shift Sf normal- 
ized by the frequency shift for a total extraction, obtained for 
single crystals #A and #C. Open circles correspond to previ- 
ous results on polycrystalline powders [l7| divided by a factor 
10. The dashed line is the T 2 law below 1.8 K reported for the 
powder. In the case of the single crystals, a better fit is ob- 
tained with an exponential law (solid line). Inset: frequency 
shift at the critical temperature for both single crystals. 



superconductors with a fully open gap : X(T) cx-^/A/fceT 
e -A/k B T^ This expression is valid for fc^T < A/5, and 
leads to A/ks = 11.6(1) K for sample #A (and B) and 
A/k B = 12.3(1) K for sample #C. Note that this fitting 
procedure can lead to a slightly overestimated A value 
(up to 10%, depending on the range of the fit and the ra- 
tio between A and T c ) but unambiguously shows that the 
gap is fully open in good agreement with previous tun- 
neling spectroscopy [2^, [2?} , NMR Q , and specific heat 
measurements which led to A/k B ~13.0(2) K, 10.5 K and 
13 K, respectively. 

However, fitting the low temperature data only leads 
to the size of the minimum superconducting gap. To 
unambiguously exclude the presence of any other gaps 
(and/or other gap symmetries) it is necessary to analyze 
the full temperature dependence of the normalized super- 
fluid density p s (T) oc 1/A(T) 2 up to T c . This superfluid 
density can be deduced : 

- either from the temperature dependence of the lower 
critical field : H cl = &l/(4n\ 2 )(Ln(K) + c(k)) where 
K = A/£ (£ being the coherence length) and c(n) a k 
dependent function tending towards ~ 0.5 for large k 
values. As n is almost temperature independent (being 
in the order of 40), H c i(T) is directly proportional to the 
superfluid density which we will call Pg cl ■ 

- or by introducing the absolute value of the pen- 
etration depth at T=0K (A ) into the TDO data : 
Ps DO {T) oc l 1+sx { T)/x J l = [ 1+f s /(r) /i/ xi?./A ] 2 where 
R is a geometrical factor|28j. 

The local magnetic induction has been measured with 
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Figure 2: Remanent field B rem in flux quantum units (<&o/S, 
S being the active area of the probe) as a function of the 
applied field H a in sample #A showing that B rem remains 
close to zero up to H a = H p (see text for details). Left inset: 
Local induction at T = 4.2 K as a function of the applied field 
for several probe positions (see right inset) showing that, even 
the penetration is much stronger close to the edges (probe 4), 
the same first penetration field (~35 G) can be obtained on 
all of the probes. Right inset: Field profiles at T — 4.2 K 
for different values of the applied field (measured on probe 2) 
clearly showing the Bean profile characteristic of bulk pinning. 
Probe 8 is located close to the center of the sample and probe 
4 close to the sample edge. The spacing between probes is 
20pm. 



a miniature 16x16 um 2 Hall probe and the first pen- 
etration field Hp has been deduced by measuring the 
remanent field (B rem ) in the sample after applying an 
external field H a and sweeping the field back to zero. 
For H a < Hp (i.e. in the Meissner state) no vortices pen- 
etrate the sample and the remanent field remains equal 
to zero (actually close to zero due to partial penetra- 
tion through the sample corners). H a is then progress- 
ively increased until a finite remanent field is obtained 
(see Fig. 2). Indeed, since vortices remain pinned in the 
sample, B rem rapidly increases for H a > H p , varying 
as (H a - H p ) a with a = 0.4 ± 0.1 (solid lines in Fig.2 
for a = 0.5). We get H p ~ 50 ± 5 G, - 55 ± 5 G and 
~ 70 ± 10 G for sample #A, #B and #C, respectively. 
In samples with rectangular cross sections, H p is then 
related to H c \ through H c \ « H v /tanh { \/ ad /2w) where 
a varies from 0.36 in strips to 0.67 in disks |23.l30l|. Tak- 
ing an average a value ~ 0.5 we hence get -ff^O) ~ 
125 ± 15 G and correspondingly Ao = 230 ± 15 nm (intro- 
ducing H c2 (0) = $ /2<(0) 2 = 9.5 T 0). 

This value is in good agreement with the value deduced 
from muon spin relaxation data in ceramics 13 if. previous 
lower critical field measurements in powder [32j as well 



as the value calculated from the thermodynamic critical 
field deduced from specific heat measurements by Walte 
et al. 

Note that, as pointed out by Walte et al., this value 
is much larger than the London clean limit BCS value 
•\l(0) = c/iOp ~ 60 nm (u> p being the plasma fre- 
quency ~ 3.2 eV [lH ). In the presence of strong mass 
renormalization and/or impurity scattering effects Ao = 
A l (0) y/1 + \e-phy/l + £o/l where / is the mean free path 
and £o ~ ^f" ( v f being the bare Fermi velocity). In- 
troducing \ e ~ph ~ 1.8 and vp ~ 2.1 x 10 5 m.s _1 fl3l |. 
one obtains I ~ £ /4 ~ 10 nm hence confirm that both 
strong coupling and strong impurity diffusion are present. 
Note that this I value corresponds to a resistivity p = 
vp/eolujp ~ lOuficm i.e. slightly lower than the resid- 
ual resistivity measured in similar crystals p ~ 30 pficm 
[20l |. However, since p is expected to be in the order of 
X^PQix/S./hil + X e -ph) (dirty limit), a residual resistivity 
of 30 pficm value would thus require that \ e - p h <S 1 in 
striking contrast with reported values. 

Pg cl (solid symbols) and p^ DO (open symbols) are 
displayed in Fig. 3 for sample (squares) and #C 
(circles). The two techniques lead to strikingly differ- 
ent temperature dependence for the superfluid density. 
For an isotropic superconducting gap, the BCS superfluid 
density ps(T) reduced by thermally activated excitations 
is expected to be given by : 

where / is the Fermi Dirac distribution, E the energy 
above the Fermi energy, A(T) the value of the super- 
conducting gap at the temperature T. As shown in 
Fig. 3 (solid lines) very good fits to the p^ cl data are 
obtained using an alpha model in which the temperature 
dependence of the superconducting gap (normalized to 
its T=0K value) has been assumed to be equal to the 
reduced BCS weak coupling value calculated from the 
gap equation [33[ and taking A(0) = 2fcsT" c . Note that 
a superconducting gap equal to its weak coupling theory 
value (A(0) = 1.76fcsT c ) only leads to a poor fit of the 
data, confirming the large value of the A(0)/fcsT c ratio 
previously obtained by bulk probes such as specific heat 
measurements. 

On the other hand, p T g DO displays a strong downward 
curvature at low temperature followed by a clear upward 
curvature as the superfluid density drops below 0.5 (i.e. 
for A(T) > 1.4Ao). As pointed out above one has to de- 
termine the R/Xo ratio in order to deduce from the 
df/Afo data. The R value has been calculated from the 
aspect ratio using the formula introduced by Prozorov 
|28| . The validity of this procedure has been checked 
on Pb samples. Moreover, different AC magnetic field 
orientations on the same single crystal of MgCNi3(but 
different R) show the same quantitative temperature de- 
pendence of A(T), consistently with an isotropic cubic 
system. 
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Figure 3: (Color online) Normalized superfluid density de- 
duced from H c i measurements (full symbols) and TDO meas- 
urements with Ao =230 nm (open symbols) for samples 
(blue squares) and #C (red circles). The solid lines are the 
fit for a superconducting gap A = 2fcsTc with T c =4.8, 5.6, 
6.2 and 7K (see text). Inset: influence of the Ao value used 
to deduce ps from TDO measurements. 



A possible explanation would hence be an underestim- 
ation of Ao- The influence of Ao is displayed in the in- 
set of Fig. 3. As shown taking Ao — 700 nm instead of 
230 nm leads to a temperature dependence for p T s sim- 
ilar to the one obtained for Pg cl - This value is however 
well above our error bars on Ao and would correspond to 
fioH c i(0) ~ 15 G i.e. even smaller than our first penetra- 
tion field values (~ 55 G). Note that strong bulk pinning 
could lead to an overestimation of H p , if measured in 
the center of the sample (see for instance [34[) but we 
checked that very similar H p values are obtained for sev- 
eral probe positions by placing the sample on an array of 
11 miniature (10 x 10 pm 2 ) probes : as shown in the left 
inset of Fig. 2, the field distribution clearly presents the 
V — shape profile characteristic of a strong bulk pinning. 
Even though those profiles confirm the good homogen- 
eity of the sample, one can not exclude the presence of 
a strong disorder at the surface of the samples leading 
to a surface penetration field much larger than the bulk 
value. However, the Ao=700nm value would require an 
extremely small mean free path (~ 1 nm, see discussion 
above). A possible difference between the mixed state 
and Meissner state penetration depth values associated 
either to a Doppler shift induced by the supercurrents on 



the excitation spectra [35l . |36| or to a strong field depend- 
ence of A in the mixed state (see for instance [37]]) due to 
multiband effects can be excluded in our isotropic, fully 
gapped system. 

Another explanation could be a difference between bulk 
and surface critical temperature. Indeed, at low temper- 
ature TDO measurements only probe the sample on a 
typical depth in the order of Ao ~ 0.2 urn, i.e. roughly 
0.4% of the total volume (for a Volume to Surface ratio of 
50 pm) . In the presence of a weak coupling superconduct- 
ing gap, this volume only increases to about 20% of the 
sample volume for T — > T c /2 and the bulk of the sample 
is only probed close to T c as the magnetic penetration 
depth finally diverges for T — > T c . On the other hand, 
the Hall probe has been placed close to the center of the 
sample in the HPM measurements and is hence sensitive 
to the bulk of the sample. In the case of MgCNi3, it is 
known that the critical temperature has a surprising high 
sensitivity to a very small change in the C or Ni stoech- 
iometry jj, [2(| and also surface stress (Til , [H[ . Assuming 
that the critical temperature of the surface is 20 % smal- 
ler than the bulk value, very good fit to the data could 
be obtained for p^ DO using Eq.l for T < | T c (still tak- 
ing A(0) = 2kBT c , see solid lines in Fig. 3). Note that 
a large dispersion of the T c values in powder might ex- 
plain the anomalous temperature dependence observed 
in previous A measurements. Clear deviations from the 
standard BCS theory (Eq.l) have been observed in sys- 
tems like MgB2 [39| or more recently in pnictides [40| 
but in our case those deviations in p]^ 00 are due to sur- 
face inhomogeneities (disorder and/or T c ) and our meas- 
urements emphasize the importance of coupling comple- 
mentary experimental probes in order to unambiguously 
address this issue. 



To conclude, we have shown that the temperature de- 
pendence of the magnetic penetration depth is exponen- 
tial in MgCNi3 single crystals signalling the presence of 
a fully open superconducting gap. A drastically different 
behavior has systematically been observed between the 
superfluid density extracted from the lower critical field 
and TDO measurements performed on the same sample, 
which are most probably due to surface disorder and/or 
a depletion of 20% of the critical temperature at the sur- 
face. 
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